Semileptonic decays of D mesons in three-flavor lattice QCD 
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We present the first three-flavor lattice QCD calculations for D —* nlu and D ^ Klu semileptonic decays. 
Simulations are carried out using ensembles of unquenched gauge fields generated by the MILC collaboration. 
With an improved staggered action for light quarks, we are able to simulate at light quark masses down to 
1/8 of the strange mass. Consequently, the systematic error from the chiral extrapolation is much smaller 
than in previous calculations with Wilson-type light quarks. Our results for the form factors at g'^ =0 are 
/^'""'(O) = 0.64(3)(6) and /+ ^^'(0) = 0.73(3)(7), where the first error is statistical and the second is 
systematic, added in quadrature. Combining our results with experimental branching ratios, we obtain the CKM 
matrix elements \Vcd\ = 0.239(10)(24)(20) and \Vcs\ = 0.969(39)(94)(24), where the last errors are from 
experimental uncertainties. 



PACS numbers: 11.15.Ha,12.38.Gc,13.20.Fc 

Semileptonic decays of heavy-light mesons are of great in- 
terest because they can be used to determine CKM matrix el- 
ements such as \Vub\, \Vcb\, iKdl and \Vcs\- The accuracy of 
one of the most important, \Vub\, is currently limited by large 

n 

theoretical uncertainty 1 1]. Lattice QCD provides a systemat- 
ically improvable method of calculating the relevant hadronic 
amplitudes, making the determination of | Vub \ and other CKM 
matrix elements more reliable and precise. 

Semileptonic D meson decays, such as Z) — > Klv and 
D — > Ttlv, provide a good test of lattice calculations, because 
the corresponding CKM matrix elements | Vcs \ and | Vcd \ are 
known more accurately than \Vub\ The decay rates and 
distributions are not yet very well known, but the CLEO-c ex- 
periment plans to measure them with an accuracy of a few per 
cent. Furthermore, measurements of leptonic and semilep- 
tonic D(j,) decays can be combined so that the CKM matrix 
drops out, offering a direct and stringent check of lattice QCD. 

Recently, dramatic progress has been achieved in lattice 
QCD, for a wide variety of hadronic quantities. Reference |2] 
showed agreement at the few per cent level between three- 
flavor lattice QCD and experiment for f.^, Jk, mass split- 
tings of quarkonia, and masses of heavy-light mesons. The 
main characteristics of these quantities are that they have at 
most one stable hadron in the initial and final states, and 
that the chiral extrapolation from simulated to physical light 
quark masses is under control. This class can be called "gold- 
plated" and many of the lattice calculations needed to test 



the Standard Model are in this class. The work reported here 
is part of a systematic effort to calculate the hadronic matrix 
elements needed for leptonic and semileptonic decays, and for 
neutral meson mixing |3l|4|]. 

In this Letter we report results for D Klv and D ^ nlv 
semileptonic decay amplitudes. All previous lattice calcula- 
tions of heavy-light semileptonic decays have been done in 
quenched = 0) QCD. In addition to quenching, they also 
suffered from large uncertainties from the chiral extrapolation 
and, in some cases, from large heavy-quark discretization ef- 
fects. Here we bring all three uncertainties under good-to- 
excellent control. Indeed, this Letter presents the first calcula- 
tion in unquenched three-flavor lattice QCD, where the effect 
of dynamical u, d and s quarks is correctly included. 

The relevant hadronic amplitude {P\V^\D) {P = tt, K) is 
conventionally parametrized by form factors /+ and /o as 

{P\V^\D) = f+{q^){pD+pp-^Y + Uq')^'' (1) 

where q — ~ pp, — {rnj^ — mp) q^/q^. The differ- 
ential decay rate dT/dq^ is proportional to |Vca:p|/+(g^)p, 
X — d,s. (A contribution from /o is proportional to the lepton 
mass squared.) We calculate /+ and /o as a function of q^ 
and determine the decay rate F and the CKM matrix | Vex \ by 
integrating over q'^. Preliminary results have been 

reported in Ref. |4, 5]. 

Our calculations use ensembles of unquenched gauge fields 
generated by the MILC collaboration |6] with the "Asq- 
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tad" improved staggered quark action and the Symanzik- 
improved gluon action |7]. The resuhs are obtained on 
the "coarse" ensembles with sea quark masses amf^''^ = 
0.005, 0.007, 0.01, 0.02, and 0.03. The gauge coupHng is ad- 
justed to keep the same lattice cutoff (a^^ « 1.6 GeV) and 
volume [L^ x T « (2.5 fm)"^ x 8.0 fm]. Each ensemble has 
about 400-500 configurations. For details of the gauge con- 
figurations, see Ref. |61. 

For the light valence quarks, we adopt the same staggered 
action as for the dynamical quarks. The valence light (u, d) 
quark mass to™' is always set equal to m^'^^. The valence 
strange quark mass is oto™' = 0.0415, which is slightly 
larger than the physical value anis = 0.039 (at this lat- 
tice spacing) determined from fixing the masses of the light 
pseudoscalars |6]. We have repeated the calculations with a 
strange quark mass slightly too small, and find a negligible 
difference. Since the computation of the staggered propaga- 
tor is fast, we can simulate with to/ as low as TOs/8. Con- 
sequently we are able to reduce the systematic error from the 
chiral extrapolation (to/ mud) to « 3%, as we show be- 
low. In contrast, previous calculations with Wilson-type light 
quarks simulated at to; > to=/2 and typically had 0(10%) 
errors from this source alone |2ll • 

For the valence charmed quark we use the clover action 
with the Fermilab interpretation 1 9] . The bare mass is fixed via 
the Ds kinetic mass 1 3]. The free parameters of both the action 
and the current are adjusted so that the leading heavy-quark 
discretization effects are 0{asaAQCD) and O ((oAqcd)'^), 
where Aqcd is a measure of the QCD scale. 

The hadronic matrix element {P\V^\D) is extracted from 
the 3 -point function in the D meson rest frame (pjj ~ 0) 

C^,^''{t.,ty;p) ^e^f«(Op(0)V;(y)Ol,(x)), (2) 

where p = pp, Vfj. = "ipcltiipx (x = d, s) is the heavy-light 
vector current on the lattice, and On and Op are interpo- 
lating operators for the initial and final states. The heavy- 
light bilinears Vf^ and Od are formed from staggered light 
quarks and Wilson heavy quarks as in Ref. 1 10]. The 3-point 
functions are computed for light meson momentum p up to 
27r(l, 1, 1)/L, using local sources and sinks. The sink time is 
fixed typically to = 20. To increase the statistics, the cal- 
culations are carried out not only at the source time to ~ 
but also at to = 16, 32, 48 (and and ty shifted accord- 
ingly). The results from 4 source times are averaged. Sta- 
tistical errors are estimated by the jackknife method. To ex- 
tract the transition amplitude {P\V^\D) we also need meson 
2-point functions Of (t,;p) = Ex e'^ "(OM(0)Ol,(a:)), 
where M — D,7r,K. They are computed in an analogous 
way. For the light meson (M = tt, K) the 2-point function 
couples to the Goldstone channel of staggered quarks. 

A drawback of staggered quarks is that each field pro- 
duces 4 quark species, called "tastes" to stress that the extra 
3 are unphysical. One consequence is that the light meson 
2-point function contains a 16-fold replication of the desired 



hadrons. On the other hand, the heavy-light 2-point function 
C2 does not suffer from such replication, because contribu- 
tions of heavy quarks with momentum p ^ 0(7r/a) are sup- 
pressed 1 10]. The same holds for 3-point functions that in- 
clude at least one Wilson quark, such as C^^^ . To check 
these properties, we carried out a preparatory quenched cal- 
culation ] 4 ] , finding reasonable agreement with those obtained 
previously with Wilson light quarks fs^ . 

Another consequence of the extra tastes is that the 3- and 2- 
point functions receive contributions from states that oscillate 
in time, in addition to the ground state and non-oscillating ex- 
cited state contributions. For example, the 3-point function's 
time dependence takes the form 

+ (-l)*«^ie"-^'*''e^-^"(*--*«) + • • • ,(3) 

whereto cx (P|V"'"|£i). 

As usual, the desired hadronic amplitude is extracted from 
fitting the 3- and 2-point functions. We employ two meth- 
ods. In the first method, we form the ratio R{ty) = 
C^^'^^ {tx.ty) l[C2 ['ty)C2 [tx - ty)], and fit to a constant in 
ty. The oscillating state contributions are partly canceled in 
the ratio, and further reduced by taking the average, R{ty) — 
[R{ty) + R{ty + l)]/2. A plateau is then found for ty around 
tx/2. In the second method, we first fit C^^^ and O^'^ sep- 
arately, using a multi-exponential form similar to Eq. (|3}, and 
then obtain {P\V^\D) from the fit results. The results from 
the two methods always agree within statistical errors. The 
difference between two results is less than 3% for the lower 
two momenta, and as large as 3% for the higher two momenta. 
We choose the first method for central values and take 3% as 
the systematic error from the fitting. 

The lattice heavy-light vector current must be multiplied 
by a renormalization factor Zy . We follow the method 
in Ref. writing Z{>^^ = pv^{Z^^ Z^"")^/^ . The flavor- 
conserving renormalization factors Zy and Zy-^ are computed 
nonperturbatively from standard charge normalization condi- 
tions. The remaining factor pv is expected to be close to 
unity because most of the radiative corrections are canceled in 
the ratio ]11]. A one-loop calculation gives fl?] pv^ ~ 1.01 
and ~ 0.99 which we use in the analysis below. This per- 
turbative calculation is preliminary, but it has been subjected 
to several non-trivial tests. 

Rather than calculating the conventional form factors /o 
and /+ directly, we first extract the form factors /|[ and f±, 
as in Ref. ||8||], defined through 

{P\V'^\D) = V^[v^ h{E)+p^^fdE)\, 

where v = po/mD, P± = Pp — Ev and E = v ■ pp is the 
energy of the light meson. /|| and f± are more natural quanti- 
ties in the heavy quark effective theory, and chiral expansions 
are given for them as a function of E lll3lfl4ll . We therefore 
carry out the chiral extrapolation in to/ for /|| and f± at fixed 
E, and then convert to /o and /+. 
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FIG. 1: a-'/^x 

as a function of {aE) for the D ~* tt decay. Sym- 
bols are raw data and lines are fitting curves with the parametrization 
of Eq. 0. Results at m; = 0.03, 0.02 and 0.01 are shown. 



To perform the chiral extrapolation at fixed E, we inter- 
polate and extrapolate the results for /y and f± to common 
values of E. To this end, we fit /y and f± simultaneousl y us - 
ing the parametrization of Becirevic and Kaidalov (BK) 11511 . 



(l-g2)(l_ag2)' 



(4) 



where (f — /m\,,, and F = /+(0), a and j3 are fit param- 
eters, and /+, /o and are converted to /y, /j^ and E before 
the fits. An advantage of the BK form is that it contains a pole 
in f+{q^) at q^ = m\),, where toc is the lattice mass of the 
charmed vector meson with daughter quark x. The BK fit for 
/j^ is shown in Fig. [2 using data for all available momenta p. 
Excluding the data for the highest momentum 27r(l, 1, 1)/L 
gives indistinguishable results. 

We perform the chiral extrapolation using recently obtained 
expressions 1 14] for heavy-to-light form factors in staggered 
chiral perturbation theory (S^PT) |16]. Compared with the 
continuum xPT formulae 1 13], the S^PT formulae include 6 
new parameters (4 splittings and 2 taste-violating hairpins), 
which parameterize lattice discretization effects. The new pa- 
rameters are fixed from the analysis of light pseudoscalars 
10]. As usual, low energy constants appear, such as the chi- 
ral coupling / and heavy- to-light meson coupling g. We 
take / = 130 MeV and g = 0.59, but changing these con- 
stants by 10% has negligible effect. The fit form we adopt 
("SxPT-Hlinear") is 



f^^\\{E)=A[l + 5f^^\\{E)]+Bmu 



(5) 



where A, B are fit parameters, and (5/_i.y is the SxPT correc- 
tion. To estimate the systematic eiTor here, we try a simple 
linear fit and a "SxPT-nquadratic" fit with a term Cm^ added 
to Eq. Q. A comparison of the three fits is shown in Fig.|2l 
For the D ^ tt [K] decay the linear fit gives 3% (2%) larger 
results at mi = rriud- The results from the S^PT-nquadratic fit 
typically lie between the results from the other two fits, with 
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FIG. 2: m;-dependence and chiral fits for a^^^^ ff^^ for sev- 
eral values of (aE)^. The SxPT+linear fit(solid), SxPT+quadratic 
fit(dotted) and linear fit(dashed). 
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FIG. 3: D ^ TT and D ^ K form factors. The experimental values 
are taken from Ref. L17.1 . 



larger eiTors. We therefore take 3% (2%) as the systematic 
error from the chiral extrapolation for the D ^ tt (K) decay. 

We now convert the results for f± and /y at m; = niud, to 
/+ and /o. To extend /+ and /o to functions of q^, we again 
fit to the form Eq. The results are shown in Fig.|3l with 
statistical errors only. We then obtain the decay rates F/j Vca; P 
by integrating (phase space) x over q^. Finally, we 

determine the CKM matrix elements | Vcd \ and | Vcs \ using ex- 
perimental lifetimes and branching ratios fP]. These main re- 
sults are summarized in Table U 

The results presented above rely on the dependence of 
BK parametrization, Eq. ©. To estimate the associated sys- 
tematic error, we make an alternative analysis without it. We 
perform a 2-dimensional fit in (m/ , E) to the raw data em- 
ploying a polynomial form plus the SxPT correction <5/y,j^. 
The result from this fit agrees with the one from the fit with 
Eq. (|5j within statistical eiTors. The deviation between the 
two fits is negligible at q^ ^ 'Zmax ^^'^ about laatq^ ^ for 
f±\\, giving a 2% difference for the CKM matrix elements. 
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TABLE I: Fit parameters in Eq. J4}, decay rates and CKM matrix 
elements. The first errors are statistical; the second systematic; the 
third experimental. 

P F a 13 r/IVe^l' [ps-^] 114^1 

TT 0.64(3) 0.44(4) 1.41(6) 0.154(12)(31) 0.239(10)(24)(20) 

K 0.73(3) 0.50(4) 1.31(7) 0.093(07)(18) 0.969(39)(94)(24) 



where the first errors are statistical, the second systematic, 
and the third from uncertainties in the CKM matrix elements. 
We do not assume any cancellation of errors in the ratios, al- 
though some may be expected. Our results agree with recent 
experimental results, /^'^''(O) ^ 0.73(15), /f~*-^(0) = 
0.78(5) id, /^'^''(0)//f ^-^(0) = 0.86(9) and r(L>" ^ 
TT-e+v^)lT{D^ ^ K-e+v,) = 0.082 ± 0.008 



With only one lattice spacing, the systematic error from dis- 
cretization effects can be estimated only by power counting. 
The leading discretization errors from the Asqtad action are 
O (as(aAQCD)^) ~ 2% (after removal of taste-violating ef- 
fects with SxPT), taking Aqcd = 400 MeV and as = 0.25. 
In addition, there is a momentum-dependent error from the 
final state. The BK parameters are determined by the lower 
momentum data; in particular, the fits are insensitive to the 
highest momentum 27r(l, 1, 1)/L. Therefore we estimate this 
effect to be 0(as(ap)^) w 5%, taking the second-highest mo- 
mentump = 27r(l, 1, 0)/L. The HQET theory of cutoff ef- 
fects llSl Il9ll can be used to estimate the discretization error 
from the heavy charmed quark. In this way, we estimate the 
discretization error to be 4-7%, depending on the value cho- 
sen for Aqcd (in the HQET context). This is consistent with 
the lattice spacing dependence seen in Ref. fs"]. In future work 
we expect to reduce and understand better this uncertainty, so 
we shall adopt the maximum, 7%, here. 

A summary of the systematic errors for the form factors 
/+.0 or the CKM matrix elements \Vcx\ is as follows. The 
error from time fits is 3%; from chiral fits, 3% (2%) for 
D ^ TT (K) decay; from BK parametrization, 2%. The 
1-loop correction to pv is only 1%, so 2-loop uncertainty 
is assumed to be negligible. The uncertainty for a^^ is 
about 1.2% |6|; this leads to a 1% error for \Vcx\ (but not 
for the dimensionless form factors), from integrating over 
to get r/jVcxP- Finally, we quote discretization uncertain- 
ties of 2%, 5%, and 7%, from light quarks, the final state 
energy, and the charmed quark, respectively. Adding all 
the systematic errors in quadrature, we find the total to be 
[3% + 3% (2%) + 2% + 1% + 2% + 5% + 7%] = 10%. 

Incorporating the systematic uncertainties, we obtain 



/f^^(O) = 0.64(3)(6), 
f^-^iO) = 0.73(3)(7), 



(6) 
(7) 



and the ratio /f^'^(0)//f^-^(0) = 0.87(3)(9). Our re- 
sults for the CKM matrix elements (Table U are consistent 
with Particle Data Group averages |Vc<j| = 0.224(12) and 
\Vcs\ 0.996(13) fl'l; also with \Vcs\ = 0.9745(8) from 
CKM unitarity. If we instead use these CKM values as inputs, 
we obtain, for the total decay rates, 

r{D" ^ TT-l+i^) = (7.7 ± 0.6 ± 1.5 ± 0.8) x 10-3ps-\ 
r(L>" ^ R-l+v) = (9.2 ± 0.7 ± 1.8 ± 0.2) x lO^^ps"^ 
T{D^ Ti-l+v) 



This Letter presents the first three-flavor lattice calculations 
for semileptonic D decays. With an improved staggered light 
quark, we have successfully reduced the two dominant un- 
certainties of previous works, i.e., the effect of the quenched 
approximation and the error from chiral extrapolation. Our re- 
sults for the form factors, decay rates and CKM matrix, given 
in Table|l]and Eq. (|8} are in agreement with experimental re- 
sults. The total size of systematic uncertainty is 10%, which 
is dominated by the discretization errors. To reduce this er- 
ror, calculations at finer lattice spacings and with more highly- 
improved heavy-quark actions are necessary; these are under- 
way. Finally, unquenched calculations of B decays such as 
B Tilv and B Dlv are in progress, and will be pre- 
sented in a separate paper. 
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V{D^ K-l+v) 



= 0.084 ± 0.007 ± 0.017 ± 0.009, (8) 
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